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ABSTRACT 



Recent observation of microlensing technique reveals two giant planets at 2.3 AU and 4.6 AU 
around the star OGLE-06-109L. The eccentricity of the outer planet (e c ) is estimated to be O.H^q'q^, 
comparable to that of Saturn (0.01-0.09). The similarities between the OGLE-06-109L system and 
the solar system indicate that they may have passed through similar histories during their formation 
stage. In this paper we investigate the dynamics and formation of the orbital architecture in the 
OGLE-06-109L system. For the present two planets with their nominal locations, the secular 
motions are stable as long as their eccentricities (e&, e c ) fulfill + e\ < 0.3 2 . Earth-size bodies 
might be formed and are stable in the habitable zone (0.25AU-0.36AU) of the system. Three 
possible scenarios may be accounted for formation of e& and e c : (i) convergent migration of two 
planets and the 3:1 MMR trapping; (ii) planetary scattering; (iii) divergent migration and the 3:1 
MMR crossing. As we showed that the probability for the two giant planets in 3:1 MMR is low 
(~ 3%), scenario (i) is less likely. According to models (ii) and (iii), the final eccentricity of inner 
planet (e;,) may oscillate between [0-0.06], comparable to that of Jupiter (0.03-0.06). An inspection 
of et, e c 's secular motion may be helpful to understand which model is really responsible for the 
eccentricity formation. 
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1. Introduction 

To date, more than 370 exoplanets are detected 
mainly by Doppler radial velocity measurements (0), 
among them 9 planets are revealed by gravitational 
microlensing (Udalski et al. 2005; Beaulieu et al. 
2006; Bennett et al. 2006, 2008; Gould et al. 2006; 
Gaudi et al. 2008; Dong et al. 2009; Janczak et al. 
2009). The use of microlensing technique to planet- 
searching is based on the idea of general relativity 
that light passing through a mass should bend as if 
it passes a lens (Mao & Paczynski 1991). The major 
advantage of the technique is that it favors to de- 
tect planets in moderate distance to the host star (~ 
3 AU), which is complimentary to the radial veloc- 
ity measurements. Among the planets detected by 
microlensing, OGLE-06-109L system is the first one 
with observed multiple planets. It is 1490 pc away 
from the Sun, with a star of ~ 0.5 M (solar mass) 
and two planets of 0.71 Mj (Jupiter mass) and 0.27 
Mj in the orbits of 2.3 AU and 4.6 AU, respectively 



(Gaudi et al. 2008). 
bital elements. 



Table 1 lists their nominal or- 
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Several features of the system make it an anal- 
ogy of the solar system, (i) The two planets have 
positions and masses similar to those of Jupiter and 
Saturn up to scale changes, (ii) The fitted eccentric- 
ity of OGLE-06-109L c is modest (0.11±^), com- 
parable with those of Jupiter (0.03-0.06) and Saturn 
(0.01-0.09) during their secular evolution, although 
the eccentricity of OGLE-06-109L b is unknown, (iii) 
While Jupiter and Saturn's orbits are closing to 5 : 2 
mean motion resonance (MMR), OGLE-06-109L b 
and c may be close to 3 : 1 MMR. (iv) All these 
four giant planets are located well outside the snow 
lines of their systems (~ 0.68 AU for OGLE-06-109L 
and ~ 2.7 AU for solar system), indicating that, un- 
like the most observed extra-solar systems with hot 
Jupiters, the migration of these four giant planets 
were not so efficient, (v) The habitable zone of the 
OGLE-06-109L system, [0.25AU-0.36AU], may have 
stable orbits. 

Considering most of the observed exoplanets 
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have close-in orbits with an average eccentricity ~ 0.2 
(Udry & Santos 2007), the investigation of formation 
scenario in OGLE-06-109L system may bridge the 
gap between the solar system and the most observed 
exoplanet systems. Another aim of the paper is to 
predict the eccentricity of OGLE-06-109L b. Due to 
the significant uncertainties in orbital determination 
of microlensing technique (<~ 10%, S. Mao, private 
communication), the orbital parameters of the two 
detected planets in OGLE-06-109L system are poorly 
known. A precise determination by other means like 
radial velocity is still impossible due to the great dis- 
tance (~ 1490 pc) between OGLE-06-109L and the 
Sun. So the investigation of their formation and dy- 
namics is helpful to reveal their orbital parameters, 
especially their eccentricities. 

The mechanism of eccentricity excitation is not 
fully understood, especially in single planet systems. 
During the early stage of planet formation, disk- 
planet interaction tends to damp the eccentricity of 
the planet. According to the linear theory, the planet 
exerts torques on a dynamically cold disk (c < rQ, 
where c, r, are the sound speed of gas, the or- 
bital radius and angular velocity of planet motion, 
respectively) mainly at the Lindblad (LR) and coro- 
tation (CR) resonances (Goldreich & Tremaine, 1979, 
1980; Ward 1988). Only non co-orbital LRs can ex- 
cite planetary eccentricity, while CRs and co-orbital 
LRs damp the eccentricity. In the linear regime, 
torques from the latter dominate the evolution, so the 
planetary eccentricity is damped, unless the planet is 
massive enough to clear the local gas disk (Lin & 
Papaloizou 1993). Two-dimensional hydrodynamical 
simulations show that the critical mass (M cr ;t) of the 
planet above which its eccentricity will be excited un- 
der disk tide is <~ 20 Mj, and M cr ;t might be reduced 
into the range of the observed extrasolar planets at 
a very low disc viscosity (Papaloizou et al. 2001). 
During the later stage of planet formation, the deple- 
tion of the gas disk will increase the eccentricity of 
a planet through the sweeping of secular resonance 
(Nagasawa et al. 2003). 

For multiple planetary systems, there are mainly 
several scenarios that will excite the eccentricities of 
the planets: 

(i) Convergent migration and resonance trap be- 
tween two planets. A planet in a gaseous disk will 
migrate inward cither due to the imbalance of LR 
torques or co-evolute with the viscous disk when it 



is massive enough to open a gap around it. In the 
case that the migration speed of inner planet is slower 
than that of the outer one, or the inner one is stalled 
due to the clear of nearby gas, a trap into MMR be- 
tween the two planets is possible, which may result in 
the increase of both eccentricities (Lee & Peale 2002; 
Kley 2003). Due to the trap of MMR, their eccen- 
tricities remain oscillating around moderate values 
(<~ 0.1 — 0.3) at the end of evolution. The config- 
urations of GJ 876 b-c in 2:1 MMR and 55 Cnc b-c 
in 3:1 MMR are believed to be formed in this way. 

(ii) Planetary scattering. During the formation 
stage of planets, protoplanetary cores may undergo 
close encounters with the planets, causing ejections 
of the cores and the eccentricities excitation for the 
survival planets. Secular interactions between the 
survival planets (m& and m c ) may result in oscilla- 
tions of eccentricities between and a finite value 
(~ 0.1). Such a motion is near the separatrix of libra- 
tion and circulation of difference of perihelion longi- 
tude (Azubc) in the eccentricity plane (e(,e c cos Azub c , 
ei,e c smAwbc), so it is called a near-separatrix mo- 
tion (Barnes & Greenberg 2006, 2008). This model 
can account for the eccentricity properties of the v 
Andromedae system (Ford et al. 2005). 

(iii) Divergent migration and MMR crossing un- 
der interaction with planetesimal disk. After circum- 
stellar disk depletes due to disk accretion, photoevap- 
oration or planet formation within ~ 3 Myrs (Haisch 
et al. 2001), planets may undergo migration through 
angular momentum exchanges with residue embryos 
and planetcsimals (Fernandez & Ip 1984; Malhotra 
1993; Hahn & Malhotra 1999). In the solar system, 
numerical simulations show that, Jupiter will drift 
inward, while Saturn, Uranus, Neptune may migrate 
outward, resulting in a divergent migration. Dur- 
ing the migration, the cross of 2:1 MMR between 
Jupiter and Saturn excites the eccentricities of four 
giant planets (Tsiganis et al. 2005), which may re- 
sult in the formation of Trojans population of Jupiter 
and Neptune (Morbidelli et al. 2005), the later heavy 
bombardment of the terrestrial planets (Gomes et al. 
2005), and the architecture of Kuiper belt (Levison 
et al. 2008). 

(iv) Slow diffusion due to planetary secular per- 
turbation. During the final stage of planet evolution 
when planets are almost formed in well separated or- 
bits, secular perturbations between them result in a 
slow increase of stochasticity of the system. This pro- 
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cedure can be approximated as a random walk in the 
space of velocity dispersion, and the resulting eccen- 
tricities of the planets obey a Rayleigh distribution, 
which agrees with the statistics of the eccentricities 
for the observed exoplanets (Zhou et al 2007). The 
difference between this and the previous planetary 
scattering scenario is that, slow diffusion model may 
occur in a much longer time span, it is effective es- 
pecially in the later stage of planet evolution when 
the planetary orbits are well separated and there is 
no violent scattering events. 

In this paper, we investigate the eccentricity for- 
mation scenarios and dynamics of the OGLE-06-109L 
system through N-body simulations, with focuses on 
the following topics: (a) the origin of the eccentric- 
ities for the two giant planets, (b) the stability of 
the present configuration, (c) the possible existence 
of planets in the habitable zone and the outer re- 
gion. The eccentricity formation scenarios revealed 
in this paper can be extended to other multiple plan- 
etary systems. The paper is organized as follows. In 
section 2, the dynamics of the present system is in- 
vestigated, with much attention paid on the stability 
of the two giant planets system, the inner and outer 
regions. Then in section 3, we study the various sce- 
narios (i, ii, iii) as mentioned above to account for the 
eccentricity formation of the two giant planets. Con- 
clusions and discussions presented in the final section. 

2. Dynamics and Stability of Nominal 
System 

Hereafter we denote mj, and m c as OGLE-06- 
109L b and OGLE-06-109L c, and the correspond- 
ing orbital elements are detached by a subscript b 
or c, respectively. As the nominal orbital periods of 
m,b, m c are close to 3 : I MMR (with period ratio 
2.8 ± 0.7, Table I), we first check the possibility of 
the present two planets in 3:1 MMR within the ob- 
servational error. Then the dynamics and stability of 
the OGLE-06-109L system are investigated. 

2.1. 3:1 mean motion resonance? 

We assume that the masses and orbital elements 
of the two planets obey normal distributions, 



with xq the nominal value of x, and a the observa- 
tional error. The distribution of eccentricity is ob- 
tained by e = \J x\ + x\, where x\, x-i are two in- 
dependent Gaussian |T]) with xq = and different 
dispersions a\, 02, and it is the standard Rayleigh 
distribution (Zhou et al. 2007) when a\ = u-x. We 
choose values of u\ and 02 so that the most prob- 
able eccentricity of m c is 0.11, with la confidence 
interval (0.07, 0.28). Also we assume e& follows a 
Rayleigh distribution with the most probable value 
(efc). The distributions of the inclinations are the 
same with those of e/2. The remaining angles are 
randomly generated. 

We carry out 15000 runs of simulations by in- 
tegrating the full motion of three bodies (to*, to;,, 
to c ) in the three-dimensional physical space up to 0.1 
Myrs. By checking whether any one of six resonant 
angles (3\ c — Xb — i'tui,— jzu c — kfl}, — lfl c with i, j, k, I 
non negative integers and i + j + k + I = 2) librates, 
we find the probabilities that to&, to c in 3:1 MMR are 
0.82%, 2.52%, 1.88% for e b = 0.04, 0.1, 0.2, respec- 
tively. Thus the probability for mb and m c in the 3:1 
MMR is small within observational errors. 
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Fig. 1. — Various motions and orbital crossing timescale 
T c for orbits initiated from the representative plane 
(e^ cos Aro;, c , e c ). In the right (left) half-plane, the initial 
Ato(, c = (Ai37j, c = 7r, respectively). The dashed curves lo- 
cate the equilibriums of the secular system with Aroj c (t) = 
or Ara;, c (t) = n. The thin solid lines separate the domains 
that Ato5 c (4) librates about 0/n from circulation. The thick 
solid line in the upper-left corner shows the boundary above 
which 0(,(l+e(,) < a c (l — e c ). The region in red is stable (with- 
out orbital crossing) in full three-body integration up to 10 8 
years. The regions in blue are unstable with T c < 10 4 years. 
Between them are chaotic regions with T c ~ 10 4— 8 years. Mo- 
tions originating from points with notation (a-c) are plotted in 
Fig. 2. 
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Table 1: Orbital parameters of OGLE-06-109L planetary system. 



Planet 


a 


P e 


■uj m 




(AU) 


(days) 


(Mj) 


b 


2.3 (±0.2) 


1825 (±365) 


0.71 (±0.08) 


c 


4.6 (±0.5) 


5100 (±730) 0.11 (-0.04+ ' 17 ) 


0.27 (±0.03) 



from Gaudi et al. (2008) and http://exoplanet.eu 



To investigate the secular dynamics of the two 
planets, we adopt the method of representative plane 
of initial conditions (Michtchenko & Malhotra 2004). 
Due to the existence of four center of mass integrals, 
the planar three-body (the star and two planets) sys- 
tem is a Hamiltonian one with four degrees of free- 
dom, with the Hamiltonian function (e.g., Laskar & 
Robutel 1995): 



H 



i=b 



2L? 



G 



m b m c 



m b m c 



ni,. 



{x b X c 



where \i{ — G(m* + nrii), m\ — niim*/(m* + mj, 
Li = m^^/JIiOi , with a^, e^, r, being the semi- major 
axis, eccentricity, relative position vector of mb or 
m c , respectively. An averaged system is obtained by 
eliminating the short periodic terms, 



1 



(2tt)2 



Hd\\ c , 



(3) 



where At,, A c are the longitudes of mean motion of 
mi, and m c respectively. Due to the D'Alembert's 
rule, only Awbc = ~ appears in the average 
Hamiltonian, thus the secular system is integrable 
with one degree of freedom in Awbc and its conju- 
gate momentum Kb = Lb ^1 — y/l — e^j , and Av3b c 
either circulates or liberates around or tt. So the in- 
formation of secular motion with different initial con- 
ditions can be presented in the representative plane 
(eb cos Awbci e c), where Awbc is fixed at either or 

TT. 

Fig.l shows the various motions in the represen- 
tative plane of the OGLE-06-109L system, with the 
semi-major axes of the planets being the nominal val- 
ues (Table 1). Two families of the equilibriums of the 
secular system ([3]), with Awb c {t) = or 7r are plot- 
ted as dashed lines. Nearby orbits are with Awb c {t) 
librating around or tt. Between the two libration 
regions are orbits with Att7£, c -circulating, bounded by 
thin solid lines. We also plot the orbital crossing 



time (T c ) of orbits originating from the representa- 
tive plane by full three-body integrations (with other 
angles randomly chosen). T c is defined as the min- 
imum time that either cib > a c occurs or the sep- 
aration of mb and m c is smaller than one mutual 
Hill radii. From Fig.l, we see that apsidal alignment 
{Awbc ~ 0) between mj, m c can stabilize the inter- 
acting planets. 
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Fig. 2. — Some typical evolution of orbits originating from 
the representative plane of Fig.l. Panel (a) shows the motion 
with Aroj, c librating around 0. Panel (b) is the motion of 
which Ac7j, c librates around it. Panel (c) is an example of 
near-separatrix motion: is periodically back to ~ , which 
is similar to the orbits of the Jupiter-Saturn system in Panel 
(d). 

Some typical motions originating from Fig.l are 
plotted in Fig. 2. Orbits (a), (b) are from the libration 
region so that Awbc librates around or tt, respec- 
tively. Orbits (c) are with Awbc circulating and e c is 
in a near-separatrix motion, i.e., oscillates between 
and a finite value (~ 0.1) periodically, which is very 
similar to the Jupiter-Saturn system in panel (d). 



2.3. Stable regions of the two-planet system 

We integrate the three-body system with a 
second-order WHM code (Wisdom & Holman 1991) 
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Fig. 3. — Orbital crossing time T c in ej,, e c plane for the two 
giant planets in the nominal locations (Table 1). All the other 
angles roj,, U7 C , Mf,, M c are random chosen. The solid curve 
plots Hill stability criterion given by Eq. J5j. 



from the SWIFT package (Levison & Duncan 1994). 
The time step is set as the 1/12 of the period of the 
innermost orbit. Taking initial and a c with present 
nominal values, we carry out 2500 runs of integrations 
with different initial eccentricities. All the angles in- 
cluding njfc and vd c are random chosen. The orbital 
crossing time T c , defined as the minimum time that 
ab > a c or distance of mb and m c being smaller than 
one mutual Hill radius, is presented in Fig. 3. Accord- 
ing to the results, the two-planet (m;, and m c ) system 
is stable (T c > 10 8 yrs) only if 



el + el < 0.3 2 



(4) 



holds approximately. Let us compare this to another 
commonly used stability. Hill stability requires the 
ordering of the two planets remain unchanging for all 
the time, which allows the outer planet escaping to 
infinity. Topological studies of three-body systems 
give the following sufficient criterion for two planets 
being Hill stable (Marchal & Bozis 1982; Gladman 
1993), 



2m to 



L E > 1 



3 4/3 



m b m c 



2/3, 



"pair/ III* ! (mi, + m c ) 4 / 3 

(5) 

where m to t = m*+mb+m c , m pa i r = m„,mb + m*m c + 
mbm c , L and E are the total angular momentum and 
energy of the three-body system, respectively. This 
criterion is also plotted in Fig. 3. Note that the Hill 
stability criterion gives a larger region as that of T c > 
10 8 years. This implies that, the stability defined by 



orbital crossing time T c > 10 s years is stronger than 
the Hill stability, in the sense that, our stability in 
terms of orbital crossing requires the mutual distance 
of two planets being more than one mutual Hill radius 
all the time. 



2.4. Stability of fictitious planets 

We numerically integrate the orbits of a few hun- 
dred test particles in the planetary system. Such 
simulations enable us to identify regions where low- 
mass companions can have stable orbits (Rivera & 
Haghighipour 2007; Haghighipour 2008). The test 
particles are initially located in circular orbits copla- 
nar with two planets, with their mean longitudes 
randomly set. The orbital evolution time is set as 
100 Myrs. A particle is removed from the simulation 
when its stellar distance exceeds 30 AU, or it enters 
the Hill sphere of either planet. We study two cases 
that either mb and m c are in 3:1 MMR or not. In the 
non-resonance cases, nib and m c are initially located 
at the nominal elements (Table 1) with e& = 0.06 (a 
value inferred from formation scenario in next sec- 
tion). To show the dependence of stability on Au7f, c , 
we let initially vjb = 0, Mb = 0, and w c = 0, 7r/2, tt. 
M c is randomly set. In the 3:1 MMR cases, mj, m c 
are set with nominal parameters except ab = 2.11445 
AU, e& = 0.06, mb — vj c — Mb = 0, M c = tt, so 
that mb and m c are initially in 3:1 MMR, and the 
three corresponding resonance angles, 3A C — Ah — 2wb, 
3A C — Ab — vob — tu c , 3 A c — Ab — 2zu c , liberate around 
tt, 0, 7r, respectively. 
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Fig. 4. — Lifetime of test particles in OGLE-06-109L sys- 
tem at different semi-major axes. Mean motion resonances are 
marked on the top (e.g. 5:1b means the 5:1 MMR with m^). 
The grey band in (a) shows the extension of the habitable zone. 



Inner region. Fig. 4a shows the survival time of 
test particles in inner system with different initial 
semi-major axes. When nib and m c are not in 3:1 
MMR, test particles with initial a < 1.5 AU are sta- 
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Fig. 5. — The largest and smallest stellar distances of an 
Earth-mass planet with different initial semi-major axes. m\y 
and m c are initially put at the nominal locations and eccentric- 
ities (e b = 0.06) except the 3:1 MMR case, where a b = 2.11445 
AU. The grey band is the extension of the habitable zone. 



tion of the HZ. The great variation of Q and q at 
a = 0.28 AU — 0.32 AU is due to the secular res- 
onance of m b and m c (Malhotra & Minton 2008). 
Interesting to note that, the variation of Q and q 
are small when Am bc = 7r, i.e., the eccentricity ex- 
cited by the secular resonance of m bl m c is small (Mi- 
gaszewski et al. 2009). Also, when mj, m c is in 3:1 
MMR, HZ is not in their secular resonance region, so 
the variation of Q and q are also small. No matter 
in which type of Aw bc , most of orbits are in HZ ex- 
cept a G [0.28 AU, 0.32 AU] with An7 bc = 0, tt/2. 
Considering the actual HZ is wider than this, OGLE- 
06-109L is a hopeful candidate system for hosting a 
habitable terrestrial planet. Simulations in next sec- 
tion indeed show the evidence for the formation of 
super-Earth planets in its HZ. 



ble in the sense of T c > 100 Myrs, except at 5:1, 
3:1, 5:2, 2:1 MMRs with m b . The locations of these 
MMRs depend on Azu bc . However, if m b and m c 
are in 3:1 MMR, the stable region is reduce a bit to 
a < 1.4 AU. 

Outer region. Fig. 4b shows the survival time of 
test particles in the outer system. When m b and m c 
are not in 3:1 MMR, the test particles are stable as 
long as a > 9.7 AU. However, If the two planets are 
in 3:1 MMR , the stable region is enlarged to a > 7.5 
AU except the 1:3 MMR with m c . It is interesting to 
note that, the 3:1 MMR between nib and m c increases 
the stable region in the outer regions. 

Habitable zone. Kasting et al. (1993) esti- 
mated the width of the habitable zone (HZ), where 
an Earth-like planet can have liquid water on its sur- 
face, around main sequence stars. For the star with 
mass 0.5 M Q , the most conservative position of HZ is 
0.25 AU-0.36 AU, and the actual HZ could be much 
wider than this extension, e.g., the outer edge of HZ 
could be 0.36 AU-0.47 AU. From Fig.4a, if the two 
giant planets are out of 3:1 MMR, test particles with 
initial semi-major axes in HZ are stable, although 
their eccentricities will be excited by secular reso- 
nance due to m b and m c (Malhotra & Minton 2008). 
To see whether the planet can maintain its orbit in 
the HZ, we put a planet with Earth mass in different 
initial locations. Fig. 5 shows the largest and small- 
est stellar distances (Q and q ) of the Earth mass 
planets. The shaded area is the conservative estima- 



3. Formation Scenarios of Eccentricity 

In this section, we test the previously mentioned 
three scenarios for the eccentricity excitation be- 
tween the two giant planets in OGLE-06-109L sys- 
tem: (i) convergent migration and resonance trap 
model, where the eccentricities are excited by the 
trap of 3:1 MMR during the type II migrations of 
two planets; (ii) planetary scattering model: the ec- 
centricities are generated by close encounters between 
some leftover embryos and the planets; (iii) divergent 
migration and MMR crossing model, when the ec- 
centricities are excited by the crossing of either 2:1 
or 3:1 MMR during the planetesimal-driven diver- 
gent migration. To simplify the problem, we assume 
the system has already in its later stage of formation 
so that both giant planets have already formed with 
their present masses, coexisting with tens of residue 
embryos with masses in the range of 0.1 Mq — 10 M® 
(Earth mass). 

For a planet embedded in a geometrically thin 
and locally isothermal disk, angular momentum ex- 
changes between the planet and the gas disk will 
cause a net momentum lose on the planet, which 
results in a fast and so called type I migration of 
the planet dGoldreich fc Tremaine 1979HWard 1997 : 
Tanaka et al. 20021 ) . Some mechanisms are proposed 
recently to reduce the speed or even reverse the di- 
rection of migration. Laughlin et al. (2004), Nel- 
son & Papaloizou (2004) proposed that, in the lo- 
cations where the magnetorotation instability (MRI) 
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is active, gravitational torques arising from megne- 
tohydrodynamical turbulence will contribute a ran- 
dom walk component to the migratory evolution 
of the planets, thus prolong the drift timcscale. 
Paardekoopcr & Mellema (2006) noticed that the in- 
clusion of radiative transfer can cause a strong reduc- 
tion in the migration speed. Subsequently investiga- 
tions (Barutcau & Masset 2008; Kley & Crida 2008; 
Paardekoopcr & Papaloizou 2008) indicated that the 
migration process can be slowed down or even re- 
versed for sufficiently low mass planets. Through full 
3D hydrodynamical simulations of embedded planets 
in viscous, radiative discs, Kley et al. (2009) con- 
firmed that the migration can be directed outwards 
up to planet masses of about 33 M®. Due to the 
vagueness of type I migration, we do not consider this 
effect at the present paper. A detailed study of forma- 
tion of Earth-like planets in OGLE-06-109L system, 
which includes the type I migration of embryos, will 
be presented in a subsequent paper (Wang & Zhou, 
in preparation). 



3.1. Disk model 

According to the conventional core accretion sce- 
nario of planet formation, planet formed through 
planetcsimals coagulation by means of runaway 
growth and became protoplanetary embryos by oli- 
garchic growth in the protoplanetary disk (Safronov 
1969; Kokubo & Ida 1998). To model the masses of 
embryos formed in disk, we adopt the empirical min- 
imum mass solar nebula (hereafter MMSN, Hayashi 
1981) so that the surface density of gas disk at stellar 
distance a is given as 



S 9 = 2.4 x i&f gM JL.)-*/* 



g cm 



(6) 



where f g is the gas enhancement factor, /d ep = 
exp(— t/rdep) is the gas depletion factor due to disk 
accretion, photoevaporation or planet formation with 
a timescale of Td cp ~ 3 Myrs (Haisch et al 2001), t is 
the evolution time. The surface density of solid disk 



is given as, 



S d = 10/ d7icc ( 



1AU 



)- 3/2 gcm- 2 , 



(7) 



where fd is the solid enhancement factor, 7j cc is the 
volatile enhancement with a value of 4.2 or 1 for ma- 
terial exterior or interior to the snow line (0.68 AU for 
OGLE-06-109L system), respectively. In such a disk, 



the embryos will grow under cohesive collisions in a 
timescale of (Kokubo & Ida 2002; Ida & Lin 2004) 



yr, 



(8) 



where M c is the core mass. The core growth will 
continue until it accretes all the dust material round 
its feeding zone (<~ 10 Hill radii) so that an isolation 
body is achieved with mass of (Ida & Lin 2004) 



■1AU 



For the OGLE-06-109L system, the isolation mass at 
4 AU with fd = 2 is around 12 M®, above the critical 
mass (~ 10 M®) for the onset of efficient gas accretion 
to form giant planets (Pollack et al. 1996), and the 
core growth timescale is ~ 1.2 Myrs. 

Interactions between embryos and the gas disk 
may damp the eccentricities of the embryos (Gol- 
dreich & Tremaine 1980). The timescale of the 
eccentricity-damping for an embryo with mass m can 
be described as (Cresswell & Nelson 2006), 



,e Qe_ r M, w M* 

( e jcmb " 0.78 ( m 



->(,2 



a 2 S, 



.)^fn-\l + \{e r -f] yr, 



(10) 

where r, e, h, ft are the stellar distance, eccentricity 
of the embryo, scale height of the disk and the Kepler 
angular velocity, respectively, Q e = 0.1 is a normal- 
ization factor to fit with hydrodynamical simulations. 

As an embryo grows to a massive planet (> 
30 M®), it will induce strong tidal torques on the 
disk to open a gap around it (Lin & Papaloizou 1993). 
Then the planet will be embedded in the viscous disk 
to undergo type II migration. The timescale of type 
II migration for a planet with mass m p can be mod- 
elled as (Ida & Lin 2004) 



m = 



— = 0.8xl0 6 yr/"^ 



M Jup )( M Q > 



TO 



1AU 



-1/2 
f/\ (11) 



where a is a dimensionless parameter to adjust 
the effective viscosity, and we set a = 10~ 4 as a stan- 
dard value in our simulation. When a giant planet is 
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embedded in a gas disk, tidal interaction of disk may 
damp its eccentricity if it is not massive enough. As 
we mentioned in the abstract, the situation is quite 
elusive for different mass regime of the giant planets, 
so we adopt an empirical formula (Lee & Peale 2002) 



(;)- 



-K 



(12) 



to describe the eccentricity-damping rate of giant 
planets in the OGLE-06-109L system, where K is 
a positive constant with a value ranging 10 — 100 
(jSandor et al. 20071 ). After some tests, we choose 
K = 10 in this paper to let eft and e c have reasonable 
convergent values. 



3.2. Numerical simulations 

In this section, we simulate the configuration 
formation for the OGLE-06-109L system with N- 
body models. We assume the two giant planets have 
formed with the observed masses at 4 AU and 8 AU-9 
AU respectively. The physical epoch corresponds to 
this assumption is ~ 1 Myrs after the formation of 
star, so that the gas disk is still present. The two 
giant planets have opened gaps around them and will 
undergo type II migration according to equation pip 
in the viscous disk. At the initial stage of our simula- 
tion, there are some leftover embryos in inner orbits 
that have obtained their isolation masses. To mimic 
the formation of Earth- like planets in inner orbits, we 
put 18 embryos, with masses ranging from 0.17 M s to 
9 Mq derived from equation ((9J and initial locations 
from 0.25 AU to 3 AU. The mutual distances among 
the embryos are set as 10 Hill radii. An additional 
embryo with in situ isolation mass will be put be- 
tween Wf, and m c in model 2. All the planets and em- 
bryos are initially located in near-coplanar and near- 
circular orbits (e = 10~ 3 , inclination i = e/2), their 
phase angles (mean motion, longitude of perihelion, 
longitude of ascending node) are randomly chosen. 
The acceleration of the planet (embryo) with mass 
rrii is given as, 

d_ y GjM* +mj) Tj 
dt 1 r, ; 2 n 



N 

■ / ■ I J ■*■ 'l 

H"F c damp( + F m igii) 



where r^, Vj are the position and velocity vectors of 
rrij in the stellar-centric coordinates, 



cdamp 



= -2 



(vi • v i )Y l 

9 

rfT e 



(14) 



is damping acceleration, effective for all embryos and 
gas giants but with different r e in equations (|10|) and 
(fl"2|) . respectively, and the acceleration that causes 
the type II migration, 



migll 



2th 



(15) 



(13) 



is adopted for two giant planets. We numerically in- 
tegrate the evolution of equation (|13[) with a time- 
symmetric Hermit scheme (Aarseth 2003). The sim- 
ulation is performed up to 10 Myrs. As we assume 
that the gas disk depletes exponentially in a timescale 
Tdep = 1 Myrs, the gas disk almost disappears at the 
end of simulation. 

During the earlier stage when gas disk is present 
(the coming models 1-2), embryos in outer disk will 
also induce a damping of giant planets' eccentricities 
through dynamic friction, which has similar effect by 
gas disk. However, as the gas disk dominates before 
gas depletion, we did not consider the presence of 
embryos in outer disk in models 1-2, except in model 
3 where embryos in outer disk are included, after the 
gas disk depletes. 

Model 1: smooth and convergent migration. In 
this model, we vary the initial locations of two gi- 
ant planets, aj, and a c . The 18 embryos with in situ 
isolation masses in inner orbits are put so that the 
outermost one is in an orbit 3.5 Hill Radii away from 
to/,. During the evolution of the typical run Rl ( 
See Table 2 for initial parameters), and m c are 
captured into 3:1 MMR at t « 1 Myrs, with three 
resonant angles liberate around either or tt with 
amplitudes around ~ 0.4 tt (Fig. 6). Their eccentric- 
ities are excited to about 0.1 inside the resonance. 
Eccentricities of the embryos in the inner orbits are 
also excited due to secular perturbations from two gi- 
ant planets, which results in their inward migration 
in the gas disk. At the end of simulation, they merge 
into 4 planets at [0.22 AU, 0.77 AU], with masses of 
4.86 M @ , 8.76 M @ , 5.79 M e , 16.06 M ffi in the order 
of increasing semi-major axes. Noticeably, the inner 
two are in the edge of the habitable zone ([0.25 AU, 
0.36 AU]) of the system. 
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Table 2: Initial parameters for two giant planets (i-pl) and embryos (i-emb) in numerical simulations of model 1 and 2 in 
section 3.2, with the outcomes of the survival embryos (f-emb) and giant planets (f-pl). P c /Pb ls the period ratios of two giants 
at 10 Myrs, the end of our simulations. 



ID 


i-pl's 


i-emb. 's 


i-emb. 's 


i-emb. 's 


f-emb 's 


f-pl's a (AU) 




a (AU) 


No. 


masses (M@) 


a (AU) 


No. 


and {P c /P b ) 


Rl 


3.8, 8.5 


18 


[0.17, 8.94] 


[0.25, 2.8] 


4 


2.53, 5.26; 3:1 


R2 


4,8 


19 


[0.17, 9.42], 16.8 


[0.26, 3], 6.5 


1 


2.40, 5.63; 3.60 


R3 


4, 8.2 


19 


[0.17, 9.42], 16.8 


[0.26, 3], 6.5 


11 


2.83, 5.02; 2.35 



The mechanism that migration of giant planets 
triggers the merge of inner embryos in MMRs or sec- 
ular resonances had been already discussed in many 
literatures, e.g., Zhou et al. (2005), Fogg & Nelson 
(2005). However, in this case, the eccentricities of 
embryos are excited by the secular perturbations of 
outside giant planets. The configuration of 3:1 MMR 
between two planets are kept to the end of the sim- 
ulation, when they are stalled at 2.53 AU and 5.26 
AU due to the severe depletion of gas disk. We did 
8 runs in this model with at, varying in [3.5 AU, 3.8 
AU], and a c in [8.2 AU, 8.5 AU]. In all the simula- 
tions, the two giant planets are trapped in 3:1 MMR, 
with their eccentricities osculating around ~ 0.1 at 
the end of evolution as in Fig. 6c. 
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Fig. 6. — Evolution of orbits in model 1 (Rl in Table 2). 
The two giants (mi, m c ) are put initially at near-circular or- 
bits with a b = 3.8 AU and a c = 8.5 AU. 18 embryos with in 
situ isolation masses are put in inner orbits. Panel (a): Evo- 
lution of semi-major axes of the 2 planets and 18 embryos. 
The grey band shows the extension of the habitable zone. The 
green dash lines represent the evolution of planet b and the red 
solid lines show the result of planet c as the same meaning in 
panel (c). Panel (b): Evolution of semi-major axis ratio of two 
giant plants. Panel (c): Eccentricity evolution of the two giant 
planets. Panel (d, e, f): Evolutions of the resonance angles 
for the 3:1 MMR between and m c . 6\ = 3A C — Af, — 2zuf,, 
62 = 3A C — Ab — roj — zu c , S3 = 3A C — A;, — 2ro c . 
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Fig. 7. — Evolution of orbits in model 2 (R2 in Table 2). The 
two giant planets (mi,, m c ) are put initially at near-circular 
orbits with = 4 AU and a c = 8 AU. The 18 embryos in 
inner orbits are put initially as in Fig. 6. An additional em- 
bryo is put at 3.5 Hill radii inside the orbit of m c . Panel (a): 
Evolution of semi-major axes of the 2 planets and 19 embryos. 
The grey band shows the extension of the habitable zone. The 
green dash lines represent the evolution of planet b and the 
red solid lines show the result of planet c as the same meaning 
in panel (c). Panel (b): Evolution of semi-major axis ratio of 
two giant plants. Panel (c): Eccentricity evolution of the two 
giant planets. 

Model 2: Planetary scattering during migration. 
Unlike the previous model that rat and m c undergo 
smooth migration, now, besides the 18 embryos in 
inner orbits, we put an additional embryo with lo- 
cal isolation mass (denoted as m c i, slightly varies for 
different locations) between the orbits of mi, and m c . 
Fig. 7 shows the results of a typical run (R2), with 
m e i = 16.8 M®. At t w 5.2 Myrs, a close encounter 
between m c and m e i occurred, which scatters m e i 
out of the system (Fig. 7). The encounter excites e c 
up to 0.1, which in turn excites eb from to ~ 0.08, 
resulting in a configuration that e;, passing in the 
eccentricity plane (e;,e c cos An7f, c , e&e c sin Avubc) al- 
most periodically, the so called near-separatrix (of 
libration and circulation of AtUf, c ) motion (Barnes & 
Greenberg 2006, 2008). The eccentricities of embryos 
in inner orbits are also excited due to the sudden in- 
crease of eb and e c , which cause strong mergers among 
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the embryos into a planet of 14.38 M s at 0.89 AU. 

If close encounters between m e i and one of the 
giant planets occur much earlier, tidal interaction be- 
tween the planets (embryos) and the gas disk will 
eventually damp their eccentricities. In run R3, the 
scattering process occurred at t = 1.5 Myrs when the 
eccentricity-damping induced by the gas disk is still 
strong (Fig. 8). As a result, ej, and e c that excited 
during planetary scattering are damped to less than 
0.01 quickly. In this case, 11 planets are left with 
their masses from 0.21 M® to 16.81 M© at [0.15 AU, 
1.10 AU]. Among them, 2 embryos are in the habit- 
able zone. 
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Fig. 8. — Evolution of orbits in model 2 (R3 in Table 2). 
The two giants (mt,m c ) are put initially at near-circular or- 
bits with a b = 4 AU and a c = 8.2 AU. The 18 embryos in 
inner orbits are put initially as in Fig. 6. An additional em- 
bryo is put at 3.7 Hill radii inside the orbit of m c . Panel (a): 
Evolution of semi-major axes of the 2 planets and 19 embryos. 
The grey band shows the extension of the habitable zone. The 
green dash lines represent the evolution of planet b and the 
red solid lines show the result of planet c as the same meaning 
in panel (c). Panel (b): Evolution of semi-major axis ratio of 
two giant plants. Panel (c): Eccentricity evolution of the two 
giant planets. 

We did 20 runs of simulations with different ini- 
tial locations of m e \ . Among them, 6 runs have close 
encounter events during the evolution, like run R2, 
resulting in similar configurations of nib and m c as in 
Fig. 7. 5 runs have earlier close encounters so that eb 
and e c are damped in the end of simulations. The rest 
9 runs do not suffer close encounters, but have some 
milder encounter events occasionally. The resulting 
eb and e c have only slight changes up to ~ 0.01. 

In a compact configuration with ab/a c > 0.48, 
semi- major axes shift to ab/a c < 0.48 by planetary 
scattering is a possible routine that could lead to the 



trap of 3:1 MMR between nib and m c . To investi- 
gate this probability, we perform additional 900 runs 
of simulations with a simplified four-body (star-two 
giants-one embryo) model. In this model, nib and m c 
are located initially at 3.7 AU and 7.6 AU so that 
(a b /a c > 0.48) , with one embryo m e \ between their 
orbits. The initial semi-major axis of m e i is chosen at 
[5.0 AU, 5.4 AU]. The mass of m e \ slightly varies at 
different locations, in the range of 13.8 M® — 14.6 M®. 
Both F e damp and Frnigii in Eqs. (fT4"|) and (fTij)) are in- 
cluded. Numerical simulations show that, 185 runs 
out of 900 ones with close encounter events occurred 
between m el and one of the planets. The epoch for 
close encounters occurred do not show a clear cor- 
relation with the relative distance between the em- 
bryo and one of the giant planets (Fig. 9). Among the 
185 runs with scattering events, 21 runs (2.3% of 900 
runs) lead to the trap of nib and m c into 3:1 MMR, 
11 runs (1.2%) result in the trap near the boundary 
of 3:1 MMR, the rest 153 runs (17.0%) do not lead to 
the trap of 3:1 MMR. We also observe 62 runs with 
nib and m c being trapped in 8:3 MMR. Fig. 10 shows 
the final eb, e c at the end of 900 runs' simulations. 
As we can see, besides those being trapped into 8:3 
MMR, eb and e c are excited significantly only in those 
185 runs with planetary scattering, with the average 
values e b ~ 0.058 and e c - 0.085 (Fig. 10). 

Model 3: Divergent migration in the pres- 
ence of planetesimal-disk. To study the effect of 
planetesimal-disk in outside orbits after the gas disk 
is depleted, we perform simulations by including the 
embryos in outside orbits but discarding those in in- 
ner orbits, since they may be in hot orbits and have 
less affections on the outer system. We set two types 
of embryos in the outer region, those with masses of 
5 M© and of 0.2 M®. After some test simulations, 
we find that the solid disk out of 10 AU has little 
effect on the evolutions of the two giant planets at 
nominal location (See Fig. 11 for a typical run). So 
we set the outer edge of solid disk within 10 AU in 
following simulations. Two groups of simulations are 
made according to different initial ab and a c . 

Group 3a: We put initially nib and m c at ab = 3 
AU and a c = 4 AU, so m c is inside the 2:1 MMR loca- 
tion (at 4.76 AU) of m&. The separation is about 3.3 
times of their mutual Hill's radii (Rijbc), above the 
threshold (2.4 Kmc), so they are Hill stable (Glad- 
man 1993) if there is no other perturbations. We 
put 63 embryos evenly at [5.5AU, 9.5AU], including 
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Fig. 9. — Orbital crossing time of the 900 simulations with 
four-body model (star-two planets-embryo). Panel (a): Or- 
bital crossing time versus the initial position of the embryo. 
Separation between the embryo and the planet b measured by 
the planet's Hill radius is marked below the figure, while that 
with the planet c is marked above the figure. Orbits without 
close encounter are represented by green stars; other symbols 
show the runs in which close encounters happen, among them 
blue squares indicate the runs in which the two planets are 
trapped into the 3:1 MMR, black circles indicate the cases in 
which the two planets are in the boundary of the 3:1 MMR, 
red circles plot the cases that the two planets are out of the 
3:1 MMR. Panel (b): the proportions of each kinds of motions 
in the total 900 runs of simulations. 
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Fig. 10. — Final eccentricities of the two giant planets in the 
end of 900 runs of simulations. If no close encounter happens 
(indicated by blue stars), the eccentricities of the two planets 
are small (~ 0.005) except 62 runs in which the two planets are 
trapped into 8:3 MMR. In the 185 runs that close encounters 
happen (represented by red circles), eccentricities are excited 
to relatively high values. The averaged ej, and e c over the 185 
runs with planetary scattering are also plotted, with the error 
bars standing for the standard deviations. 
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Fig. 11. — Evolution of a typical run in model 3 with a 
total disk mass of 27 Mgj planetesimal disk outside in [11 AU, 
16 AU]. Panel (a) and (b): evolutions of semi- major axes of 
two giant planets. Panel (c): semi-major axis evolutions of all 
planets and embryos. Panel (d): eccentricity evolutions of two 
giant planets. 



5x5 M ffi and 58 x 0.2 Mq ones, corresponding to a 
solid disk of fd = 2 with total mass of 36.6 M e in [5 
AU, 10 AU]. 
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Fig. 12. — Evolution of two giant planets in one run of group 
3a (model iii). The blue solid lines show the epoch of 5:3 
MMR, 2:1 MMR and 5:2 MMR orbital crossing between the 
two giant planets mi, and m c . Panels (a), (b), (c), (d) show 
the evolution of semi- major axes of mj, and m c , period ratio 
(Pc/P;,), and eccentricities e(,,e c , respectively. The green dash 
lines represent the evolution of planet b and the red solid lines 
show the result of planet c in panel (a), (b) and (d). 

Fig. 12 shows the evolution of mi, and m c in a typ- 
ical run, with the innermost embryos ai n = 5.5 AU. 
Under the perturbation of outer embryos, nib {mc) 
becomes unstable and undergoes inward (outward, 
resp.) migration quite soon (Fig. 12a, b). The mi- 
gration results in the quick crossing of 5 : 3 and 2 : 1 
MMR at t w 0.05 Myrs, and 0.12 Myrs (see Fig.l2c), 
the two most strong resonances between 3:2 and 3:1 
MMRs, and the system seems to be stable until m 
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(m c ) reaches around 2.82 AU (4.90 AU, resp.), with 
a drift extension of 6% (22%, resp.). Their eccentric- 
ities are excited after the crossing of 5 : 3 and 2 : 1 
MMR, with maximum values of - 0.2-0.3 (Fig.l2d). 
Finally e& and e c oscillate in [0.03, 0.15] and [0.04, 
0.19], respectively. The 5 : 3 and 2 : 1 MMR cross- 
ings of toj and m c lead to the strong scattering of 
embryos in outer orbits (Fig. 13a), which results in 
the escape of all the embryos except one with mass 
5 M® at the orbit of 15 AU and eccentricity of 0.45. 

We did 15 runs of simulations in this group, by 
changing the initial locations of the 63 embryos so 
that din of the innermost one varies in [5.2 AU, 8 AU]. 
5 : 3 and 2 : 1 MMR crossings occurred in 3 runs, with 
eccentricities e b G [0.03, 0.24] and e c G [0.04, 0.28] 
at the end of our simulations (£=10 Myrs). We did 
not observe resonance-crossings in another 3 runs up 
to 10 Myrs' evolution, with final eccentricities e b G 
[0, 0.05] and e c G [0, 0.06]. In the rest 9 runs, m b 
and m c have strong close encounter so that m c is 
scattered out of the system, with the eccentricity of 
the only survival giant planet rnj, — 0.3. The values 
of a m corresponding to these three types of outcomes 
(MMR crossing, non MMR crossing, m c being ejected 
) do not show clear correlation, which indicates the 
chaotic states of m b and m c under the perturbation 
of outer embryos. 
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Fig. 13. — Evolution of two giant planets and embryos in 
model 3. The error bar at the end of the evolution show the 
extension of a(l — e) of each survival planets and embryos. 
The green dash lines represent the evolution of planet b and 
the red solid lines show the result of planet c. Panel (a): semi- 
major axis evolutions of one run of group 3a, the blue solid lines 
indicate the epoch of 5:3 MMR, 2:1 MMR and 5:2 MMR orbital 
crossing between the two giant planets. Panel (b): semi-major 
axis evolutions of one run of group 3b, the blue solid lines 
indicate the epoch of 3:1 MMR orbital crossing between the 
two giant planets. 
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Fig. 14. — Evolution of two giant planets in one run of group 
3b. The blue solid lines show the epoch of 3:1 MMR orbital 
crossing between the two giant planets mj and m c . Panels 
(a), (b), (c), (d) show the evolution of semimajor axes of 
and m c , period ratio {P c /Pb)> an d eccentricities e;,, e c , respec- 
tively. The green dash lines represent the evolution of planet 
b and the red solid lines show the result of planet c in panel 
(a), (b) and (d). 

solid disk of fd = 2 with total mass of 35.2 M s in [5.5 
AU, 10.5 AU]. Fig. 14 shows the evolution of TO& and 
m c in a typical run, with the innermost embryos a lrl = 
6.4 AU. TOfc, m c cross 3 : 1 MMR during the divergent 
migration at t w 0.06 Myrs and t ~ 2 Myrs, causing 
the increase of e& and e c , At the end of simulation, 
m;, and m c locate at 2.27 AU and 4.71 AU, with e b G 
[0.01 - 0.07] and e c G [0.03 - 0.1]. One 5 M e embryo 
(at 11.4 AU with e = 0.024) and 4 x 0.2 M ffi ones 
(at [16 AU, 22 AU] and e G [0.08, 0.53]) are left 
(Fig. 13b). 

We did totally 41 runs (13 runs with fd — 2 and 
28 runs with lower fd) in this group. For the 13 runs 
of fd = 2, 5 out of these 8 runs with a ln < 7.8 AU are 
observed having m b , to c 's 3:1 MMR crossing, with 
final eccentricities e b G [0, 0.08], e c G [0, 0.1]. The 
rest 5 runs with a in > 7.8 AU do not have 3:1 MMR 
crossing, with final eccentricities e b G [0, 0.04], e c G 
[0, 0.04]. For the 28 runs with lower solid disks (fd = 
1, 1.2, 1.5, 1.6, 1.8) and different a in , no 3:1 MMR 
crossing is observed in runs with fd < 1.6, maybe 
due to the small mass of solid disk (M < 22 Mg). 
While for f d = 1.6, 1.8, the probability of 3:1 MMR 
crossing is — 33%. The extension of e b and e c are 
similar to that of fd = 2. 



Group 3b: m b , m c are put initially at a& = 2.3 
AU and a c — 4.6 AU, the nominal locations. We put 
56 embryos evenly at [6.4 AU, 10.4 AU], including 
5x5 M s and 51 x 0.2 M ffi ones, corresponding to a 



4. Conclusions and Discussions 

In the paper we investigate the dynamics and for- 
mation scenario for OGLE-06-109L system in terms 
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of the observed planets (nib, m c ), aiming to under- 
stand its formation history and to predict ef, that is 
not revealed by observation. According to the inves- 
tigation, we find the evolution history of nib and m c 
depends strongly on the initial conditions, i.e., the 
cores of nib and m c before efficient gas-accretion be- 
gins. 

According to the conventional core- accretion sce- 
nario of planet formation, a giant planet forms from a 
massive embryo (> 10 M©) through accreting nearby 
gas. Embryos beyond the snow line tend to have 
larger isolated masses, thus they are the ideal can- 
didates for planetary cores. However, there will be 
more than one embryo beyond the snow line. For 
example, assuming a solid disk of 2 times of the min- 
imum mass solar nebula (MMSN) for the OGLE-06- 
109L system, the space between 3 AU and 8 AU can 
be occupied by 4-5 embryos with isolation masses 
above 10 M© and mutual separations <~ 10 Hill 
radii. Due to the long quasi-hydrostatic sedimenta- 
tion stage of gas (~ several Myrs, Pollack et al. 1996), 
and the perturbation from first generation giant plan- 
ets to nearby embryos, all isolation masses may have 
the chance to grow up into second generation giant 
planets, thus the initial locations of formed planets 
can not be well determined. 

For the two giant planets in OGLE-06-109L sys- 
tem, if they formed from embryos with relatively far 
mutual distances, their initial configuration is loose, 
e.g., a\yja c < 0.48, so that they are beyond the 3:1 
MMR. Subsequent smooth migration under disk tide 
will result in a 3:1 MMR, provided suitable gas de- 
pletion timescale. We did 8 runs in model 1, all the 
simulations result in the two giant planets trapped in 
3:1 MMR, with their eccentricities being excited and 
osculating around <~ 0.1 at the end of evolutions as 
in Fig. 6c. 

If the two planets formed from embryos with a 
relatively small distance, they may have a compact 
configuration initially with at,/a c > 0.48. Then plan- 
etary scattering among residue embryos (m e i) and 
planets arc most probably the major cause of ef, and 
e c . Among the 900 simulations of model 2 we did, 185 
runs (20.5%) with close encounter events occurred be- 
tween m e i and one of the planets, with the average 
values e b ~ 0.058 and e c ~ 0.085 (Fig. 10). Only 3.5% 
of the 900 runs lead to the trap (or in the boundary) 
of mb and m c into 3:1 MMR. 



After the gas disk is almost depleted, divergent 
migration of mi, and m c caused by the residue em- 
bryos and planetesimals in outer disk may drive nib 
and m c passing through lower order MMRs. Accord- 
ing to our simulations of model 3, the crossing of 2:1 
MMR is unlikely in the OGLE-06-109L system, as it 
will excite eccentricities of nib and m c up to 0.2 — 0.3, 
and it is easy to eject m c out of the system. On 
the other hand, the crossing of 3:1 MMR is likely, 
which will excite the eccentricities up to e\> ~ 0.06 
and e c <~ 0.10. However, from our simulations, the 
required migration depends on the mass and radial 
location of the planetesimal disk. A solid disk with 
mass enhancement factor fd > 2 over the minimum 
solar nebular may be needed, with their inner edge 
within 7.8AU. Considering the similarities of solar 
system and the OGLE-06-109L system, fd <~ 2 is still 
possible for the OGLE-06-109L system with a stellar 
mass ~ 0.5 M©. 

In summary, all the three models (i) smooth, 
convergent migration and the trap of 3:1 MMR; (ii) 
planetary scattering; (iii) divergent migration and 
the crossing of 3:1 MMR, e& and e c can be excited. 
However, the probabilities, the conditions and the 
final outcomes of these three models are different. 
Smooth and convergent migration in model (i), if it 
occurs as predicted by the standard model, could re- 
sult steadily in the trap of 3:1 MMR between nib and 
m c , with eti, e c ~ 0.1 all the time. For model (ii), the 
probability of planetary scattering occurs is <~ 20%, 
they result in average e& <~ 0.06 and e c <~ 0.09, and 
nib is more likely to undergo a near-separatrix motion 
in (et,e c cos Azubc, ei,e c sinTZ7f, c ), i.e., e& passing at a 
secular timescale (~ 0.08 Myrs as in Fig. 15a). The 
probability for nib, m c 3:1 MMR crossing in model 
(iii) depends on the mass and extension of residue 
solid disk, and will result in et e [0.01 — 0.07] and 
e c G [0.03 — 0.1], but the variations in these ranges 
are in a relative shorter timescale, e.g., ~ 0.01 Myrs 
in Fig. 15b. 

Some analytical estimations of related timescale 
is helpful to reveal the different procedures corre- 
sponding to eccentricity evolution in model (ii) and 
(iii). The timescale for the secular evolution of two 
planets is given as 27r/|gi — 52 1 7 where g\, gi are the 
two cigcnfrcquencies (Murray & Dermott 1999, Zhou 
& Sun 2003). This gives 0.076 Myrs and 0.037 Myrs 
for orbits in Fig. 15a and Fig. 15b, respectively. In the 
circular restricted three-body (CRTB) framework, 
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time (Myr) 

Fig. 15. — Eccentricity evolutions of a typical run in model 
(ii) and model (iii). Panel (a): zoomed in evolution of figure 7c 
from 6 Myrs to 7 Myrs, with a b 2.40 AU and a c 5.65 AU. 
Panel (b): zoomed in evolution of figure f4d from 2.4 Myrs to 
2.6 Myrs, with a b » 2.28 AU and a c ss 4.73 AU. 



the timescale of a massless body in the 3:1 MMR of a 
perturber is estimated as 2n / (ne-\/3/i' afd(ct)) , where 
// = m'/m* is the mass ratio of the perturber, n, e 
is the mean motion and eccentricity of the massless 
body, fd(ce) is the function of Laplace coefficients of 
semi- major axis ratio a (Murray & Dermott 1999). 
Assuming uib is the perturber, e = 0.04 gives the e- 
evolution timescale of m c in mj's 3:1 MMR as 0.01 
Myrs, although the CRTB model is not a good model 
here. So the eccentricity evolution of Fig. 15a is due 
to the secular dynamics, while that in Fig. 15b is due 
to the 3:1 MMR crossing. We also find such a 3:1 
MMR timescale is kept for these two orbits up to the 
end of simulation. 

So to understand scenarios of eccentricity for- 
mation for the OGLE-06-109L system, we need more 
detailed information of their orbits. If mb and m c 
are shown by observation that they are in 3:1 MMR, 
then model (i) should be the most possible scenario. 
Based on our simulations in section 2.1, we think the 
possibility that nib, m c in 3:1 MMR is small, thus 
model (i) is unlikely. However, either model (ii) or 
(iii) can not be decided by the present observations. 
In both models, the averaged e c ~ 0.09—0.10, roughly 
agrees with the observed value 0.11^0 oL an d P re ~ 
diet the most possible value of e& » 0.06. If nib and 
m c are observed to undergo a near separatrix motion 
in (e&e c cos Awb c , ebe c smwb c ) plane in a timescale 
of secular motion, then model (ii) is favored for the 
origin of eccentricities. And we can infer that, after 
OGLE-06-109L b and c formed, either smooth migra- 



tion history during the presence of gas disk is short, 
or during their migration history, embryos may suffer 
close encounter with the planets, exciting their eccen- 
tricities. However, if e?,, e c oscillates in a timescale 
of nearby (crossed) MMR, then most probably model 
(iii) accounts for the origin of their eccentricities, and 
based on this it is possible to predict the extensions 
of residue disk mass and location through a more de- 
tailed study. 

For the stability of the OGLE-06-109L, the two 
giant planets will be stable provided e\ + e\ < 0.3 2 . 
According to the formation scenario, super-Earth 
planets may be formed inside or outside their orbits. 
Numerical simulations show the region a < 1.5 AU 
(including the habitable zone) or a > 9.7 AU is sta- 
ble. Although the habitable zone contains secular 
resonance of the system (Migaszewski et al. 2009; 
Malhotra & Minton 2008), our investigation shows 
that it is wide enough for an Earth-mass planet being 
formed and stable at least 10 Myrs. In the rare cases 
when the two giant planets are in the 3:1 MMR, the 
stable region in inner orbits is reduced to a < 1.4 AU, 
while that in the outer region is enlarged to a < 7.5 
AU. 

When extending the analyses to other multiple- 
planet systems, as close encounters between residue 
embryos are common, we expect planetary scattering 
and the consequent near-separatrix motion of eccen- 
tricities among multiple planetary systems are also 
common, which agrees with the statistics of the multi- 
ple exoplanet systems observed (Barnes & Greenberg 
2006, 2008). Divergent migration is also possible to 
sculpt the architecture of the multiple planetary sys- 
tems. Both mechanisms can account for the presence 
of modest eccentricities in the multiple planetary sys- 
tems without necessarily being trapped in MMRs. 
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10833001, 10778603), National Basic Research Pro- 
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